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A Hartree approximation is used to study the interplay of two kinds of scaling 
which arise in high-temperature superconductors, namely critical-point scaling and 



that due to the confinement of electron pairs to their lowest Landau level in the 
presence of an applied magnetic field. In the neighbourhood of the zero-field critical 

-6 

point, thermodynamic functions scale with the scaling variable (T — T c2 (B))/ B l l 2v 

O 

which differs from the variable (T — T c (0)) / B 1 ^ 2u suggested by the gaussian ap- 
proximation. Lowest-Landau-level (LLL) scaling occurs in a region of high field 
surrounding the upper critical field line but not in the vicinity of the zero-field tran- 
sition. For YBaCuO in particular, a field of at least 10 T is needed to observe LLL 
scaling. These results are consistent with a range of recent experimental measure- 
ments of the magnetization, transport properties and, especially, the specific heat 
of high-T c materials. 
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1. Introduction 

The superconducting transition in conventional low-T c materials is well described 
by the Ginzburg-Landau mean field theory. Principally because of the large cor- 
relation volume in these materials, the region in which critical fluctuations might 
be important is too small to be accessible experimentally. In high-temperature su- 
perconductors, by contrast, the critical region may be much larger. Widely varying 
theoretical estimates of the size of this region are obtained according to details of 
the criterion employed 1 , but marked deviations from mean-field behaviour have been 
observed 2-6 over a temperature range of the order of 10 K above and below T c . 

Theoretical expectations of the kind of critical behaviour which might be ob- 
served are somewhat confused. If fluctuations in the magnetic vector potential can 
be ignored, then the zero-field transition ought to be a critical point in the uni- 
versality class of the 3-dimensional XY model, as is the superfluid transition in 
4 He. When magnetic fluctuations are included, a renormalizat ion-group analysis by 
Halperin, Lubensky and Ma 7 for a (4 — e)-dimensional system reveals a runaway of 
renormalization-group trajectories, which these authors interpreted as a signal of a 
weakly first-order transition. This interpretation is confirmed by an explicit con- 
struction of the free energy 8 . On the other hand, a renormalization-group analysis 
in (2 + e) dimensions 9 indicates a second-order transition in the universality class of 
the CP N ^ model in the limit N — > 1, while a lattice simulation of Dasgupta and 
Halperin 10 is consistent with inverted XY critical behaviour. In high-T c cuprates, the 
region in which a first-order transition or inverted XY behaviour might be detected 
is probably extremely small. These are strongly type-II materials with penetration 
depths in excess of 1000 A. In this situation, renormalization-group trajectories pass 
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very close to the ordinary XY fixed point, suggesting that magnetic fluctuations can 
indeed be ignored except in a very narrow range of temperatures near T c . 

Experience of critical phenomena in, for example, fluids and magnets suggests 
that, in the presence of an applied magnetic field B, there should be a critical region 
in which thermodynamic quantities assume the scaling form A(T,B) = B aA A(x), 
where is a critical exponent associated with the quantity A and the scaling 
variable is an appropriate ratio of scaling fields. In principle, the correct scaling 
fields would emerge from a renormalization-group analysis, but this analysis is very 
difficult in the presence of an applied field. An early calculation of Prange 11 using 
the gaussian approximation suggests that B occurs in the combination B£ 2 , where 
£ is the zero-field coherence length, so that the scaling variable should be x = 
(T - T C )IB X I 2 \ where v is the coherence-length exponent. We shall argue that 
this is not quite correct, however, and that the scaling variable should be x = 
(T - T c2 (B))/B 1 / 2l/ , where the line T = T c2 (B) is a renormalized version of the 
line usually denoted by H c2 (T), the upper critical field in mean field theory. In the 
gaussian approximation, with v — -, these two scaling variables differ only by an 
additive constant. 

A different scaling form for thermodynamic functions (which is not directly asso- 
ciated with a phase transition) arises in the lowest Landau level (LLL) approxima- 
tion 12 ~ 14 , which is usually thought to be valid in the neighbourhood of the H c2 (T) 
line. Here, the appropriate scaling variable is y = (T — T c2 (B))/B l/) , where the expo- 
nent has the value = | in 3 dimensions or = | in 2 dimensions. This scaling 
behaviour is well verified experimentally for conventional superconductors 12-16 , but 
the type of scaling which applies to high-T c materials is at present a matter of some 
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controversy. It has been claimed by Welp et al that magnetization, conductivity 
and specific heat data for YBaCuO are consistent with LLL scaling. On the other 
hand, Inderhees et al 4 and Salamon et al 5 claim that experimental data is more 
nearly consistent with critical-point scaling. In fact, the predicted scaling forms 
for the magnetization and conductivity are probably too similar to be distinguished 
experimentally. For the specific heat, the LLL scaling fit exhibited by Welp et al is 
rather poor, and is achieved only by the introduction of a prefactor which has no 
theoretical basis. 

In this work, we use a Hartree approximation to study the interplay of critical- 
point and LLL scaling in an isotropic, <i-dimensional system. Although materials 
such as YBaCuO are anisotropic, layered systems, we expect that the results of 
this study should provide a reasonable qualitative guide to the scaling behaviour 
in the vicinity of the critical point, where the coherence length is much larger than 
the interlayer spacing. The Hartree approximation is in any case too crude to give 
accurate numerical estimates either of critical exponents or of scaling functions. We 
find that critical scaling may be expected in a region of small fields and temperatures 
near T c (0), while LLL scaling occurs in a region which surrounds the H c2 (T) line, 
but stops short of the zero field critical point at a field value which we estimate at 
between 10 and 100 T. The Hartree approximation is described in section 2 below, 
and a criterion for the validity of the LLL approximation is obtained in section 3. 
Sections 4-6 discuss scaling behaviour of the field-dependent coherence length, the 
specific heat and the electrical conductivity. A comparison with recent experimental 
measurements is made in section 7, and our conclusions are summarized in section 
8. 
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2. Hartree approximation 



The Ginzburg-Landau- Wilson reduced Hamiltonian may be written in a standard 
form as 



n= d d x 



|(V-*eA)0| 2 + t o |0| 2 + iA|0| 4 



(2-1) 



where, in the critical region, t can be taken as linear in temperature and A as 
a constant. We assume throughout that the magnetic flux density B = V x A 
is uniform, and equal to the applied field. Thus, the vector potential A is not a 
fluctuating variable, and the expectation value of a quantity f((f>) is 

= J Zty/(0)exp(-W)/ J V^eM-n) • (2.2) 

For a d-dimensional system, a convenient gauge choice is A = \B{— y, x,0), where 
denotes the components in the (d - 2) dimensions "parallel" to the magnetic field. 
We implement an approximation of the Hartree type by introducing an approximate 
Hamiltonian 

Ho = J d d x [\ (V - teA)(j)\ 2 + (t + fi) |0| 2 ] , (2.3) 

where ji is determined self-consistently by requiring that H and Ho have the same 
expectation value in the ensemble of Ho: 

(^-^o)o = ^(|0| 4 )o-/i(|0| 2 )o = O. (2.4) 

The expectation values are easily expressed in terms of a sum over Landau levels as 
(|0| 2 )o = /and (|0| 4 ) o = 2/ 2 , where 

_ 2eB r d d ~ 2 k 1 

4tt V ' (2vr) d - 2 k 2 + 2eBn + t + pi + eB ' 1 ' 1 

On defining 

t = t + eB + n , (2.6) 



we obtain the constraint equation (2.4) in the form 

r, A „^ r d d ~ 2 k 1 , „. 

t = t + eB + — 2eB V / - — — -— =. 2.7 

° 4tt ^ i 2tt ) d ~ 2 k 2 + 2eBn + t V ; 



The sum and integral in this expression is divergent, unless the integration is 
restricted by an upper cutoff, of the order of an inverse lattice spacing. However, 
this divergence can be eliminated by an additive renormalization of the temperature. 
On carrying out the sum and angular integrations, we find 

t = t + eB + X(2eB) d/2 - 1 f(t/2eB) , (2.8) 

where A = A/ (4ir) d / 2 , the renormalized temperature variable is t = t + \ J °° dxx~ d / 2 
and 



re 

m = / 

Jo 



dxx^ 2 



- 1 



(2.9) 



1-e" 

For 2 < d < 4, this integral is finite. We are not able to evaluate it analytically, 
but its limiting behaviour for small and large values of z — t/2eB can be obtained 
straightforwardly. When z is small, we have 

f(z) « foz d/2 ~ 2 , (2.10) 

where fo = J °° dx x 1 ~ d ^ 2 e~ x (— for d — 3). In this limit, the constraint equation 
reads 

t w t + eB + Xf (2eB)t d/2 - 2 . (2.11) 

This limit corresponds to the lowest Landau level (LLL) approximation, where the 
sum in (2.7) is approximated by the term n — 0. In this approximation, t can be 
expressed in a scaling form, which is inherited by various thermodynamic functions 12 , 
namely 

t ^ (XfoeBfr LLL (y) , (2.12) 
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where = 2/(6 — d) (so = 2/3 in 3 dimensions), y — (t + eB) / (XfoeB)^ and the 
scaling function r LLL is the solution of 

TLLL(y) = y + 2T LLL (y) 1 - 1 ^ . (2.13) 

In the opposite limit, z — > oo, corresponding to very low fields, we have 

/(z) « -/oo^ 2 - 1 , (2.14) 

where /«, = / °° &r d / 2 (l - e"*) (= 2y/n for d = 3). 

3. Validity of the lowest Landau level approximation 

One can, of course, make the approximation of neglecting all Landau levels except 
the lowest without also invoking the Hartree approximation, and perturbative calcu- 
lations based on such an approximation have been pursued by several authors 14 ' 15 ' 18 . 
One finds that the scaling form (2.12) persists, though the scaling function is no 
longer that given by solving (2.13). It is clearly essential to know where in the 
phase diagram the LLL approximation is likely to be valid, and we address this 
question within the Hartree approximation in the following way. We suppose that 
the function f(z) is well approximated by (2.10) whenever z is smaller than some 
fiducial value, say e. For z = t/2eB = e, the constraint (2.11) reads 

t = -eB[l - 2e + 2Xf (2eBe)~^] , (3.1) 

where ip = — 1=2 — d/2 (ip = 1/2 in three dimensions) and defines a locus in 
the (t,B) plane, which is shown schematically in figure 1. Points for which z < e lie 
below this line and this, therefore, is the region in which we might expect the LLL 
approximation to be valid. However, our analysis is based on the assumption that 
the order parameter (0) is negligibly small. This assumption presumably becomes 
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invalid at some distance below the mean- field H C 2(T) line t = —eB, but we are 
unable to determine whether some form of the LLL approximation survives with 
((f)) 7^ 0. Supposing that the order parameter is indeed negligible above the line 
labelled (0) ~ in figure 1 (whose location we cannot determine precisely), we are 
led to the conservative conclusion that the LLL approximation should be valid in the 
roughly wedge-shaped region labelled LLL in the figure. This region encloses a high- 
field, low-temperature portion of the line t = —eB, but stops short of the critical 
point t = B = 0. As the condition for the accuracy of the LLL approximation is 
made more stringent, by reducing the value of e, the wedge recedes to higher fields 
and lower temperatures. 

4. Scaling of the coherence length 

From the appearance of t in the propagator in (2.7), we can identify this quantity 
in terms of a temperature- and field-dependent coherence length £(t,B) as t = 
£(t,B)~ 2 . We find from the constraint equation (2.8) that it can be expressed in a 
two-parameter scaling form as 

t = (2eB)T B (6,5) , (4.1) 

where the scaling variables are 

6 = \-\t + eB){2eB)' l ' 2v , (4.2) 

5 = \-\2eBY' 2 , (4.3) 

with exponents v = l/(d — 2) and u = 4 — d, and the scaling function tb is the 
solution of 

5r B = 6 + f{r B ) . (4.4) 
8 



It is straightforward to show that r B has a double power series expansion in 9 and 
5. For small fields and temperatures close to T c , we can take the limit 5^0 with 
9 fixed, to obtain the one-parameter scaling form 

t « (2eB)r B (9, 0) = (2eB)f- 1 (-9) . (4.5) 

This should be valid in a region near the critical point t = B = 0, indicated schemat- 
ically by the shaded "Critical" region in figure 1. According to (2.10) and (2.14), 
the critical-point scaling function t b (9, 0) = f~ 1 (—9) has the limiting forms 

/ g \2/(d-2) 

t b (9,0)^1 — \ , 9^+oc, (4.6) 

/ a \ -2/(4-d) 

TB ^°)~ (-^) ' 0--oo. (4.7) 

In general, the regions where the LLL scaling form (2.12) and the critical scaling 
form (4.5) hold must be distinct. There may, however, be a crossover region in which 
both forms are approximately valid. This requires 

T B (9,0)~9- a (4.8) 

and 

tllM ~ y~ u , (4.9) 

where a = 2z/(l — 0)/(2z/0 — 1). In the Hartree approximation, we have <r = 2/(4 — d) 
and we find from (2.13) and (4.7) that there is such a crossover region, namely the 
region of very low field below T c , where 9 and y are both large and negative. It 
is, however, in this region that the approximation (0) is likely to fail, so (4.8) 
and (4.9) are not necessarily meaningful in this region. Moreover, we do not know 
whether these limiting forms of the scaling functions are valid beyond the Hartree 
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approximation, so it is not clear whether the crossover region would be accessible in 
real materials. 

In the usual way, the scaling relation (4.1) can be reformulated as 

t=\t + eB\*'T±(P,ri), (4-10) 

where (3 = 2eB\t + eB\~ 2v , rj — \t + eB\ wv and the + and — branches of the scaling 
function refer to t > —eB and t < —eB. It is straightforward to show that rf~ has a 
double power series expansion in (5 and r\. The form (4.1) is, however generally the 
more convenient. 

5. Scaling of the specific heat 

Identification of the specific heat within the Hartree approximation is somewhat 
ambiguous. On the one hand, we can define a Hartree approximation to the free 
energy density, F = — V^ 1 In [/ P0exp(— Ho)], where Ho is the approximate Hamil- 
tonian introduced in (2.3), and identify C = d 2 F /dt 2 . On the other hand, we can 
identify the entropy density of the original Ginzburg-Landau- Wilson model (2.1) as 
S = — 1(|0| 2 ) and the specific heat as C = dS/dt, and evaluate this quantity in the 
Hartree approximation: 

C = -ij>IV (5.1) 

These two definitions are not equivalent. The latter definition was adopted, for 
example, by Bray 12 and is the one we use here. It is slightly simpler, and agrees 
with the natural definition of the specific heat in the many-component limit 19 , which 
is largely equivalent to the Hartree approximation. Up to a non-universal constant 
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prefactor, the singular part of the specific heat is then given by 



C = l-|j. (5.2) 



Scaling behaviour of the specific heat now follows directly from that of t. In the 
lowest Landau level approximation, we have 

C~l-T' LLL (y) , (5.3) 

with the scaling variable y and scaling function t lll defined as in (2.12) - (2.13), 
and this naturally reproduces Bray's result 12 . Corresponding to the two-parameter 
scaling form (4.1) for t, we find 

C = 1 - \-\2eB)- a/2u C(9, 5) , (5.4) 

where a = 2 — dv = — (4 — d)/(d — 2) and 

C(9,5) = ^r B (9,5). (5.5) 

This relation between the scaling functions for the specific heat and the coherence 
length is a special feature of the Hartree approximation, as are the associated rela- 
tions 

a = 1 — 2v = —ujv (5.6) 

between the specific heat exponent a, the coherence length exponent v and the 
correction to scaling exponent u>. It seems plausible, however, that the scaling form 
(5.4) should be more generally valid. That is, we expect that a region should exist 
in which both the asymptotic critical behaviour of the specific heat and the leading 
corrections are described by a function of the form 

C = C,- C 2 B~^C (^, uB^) . (5.7) 
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In this expression, C\ and C 2 are non-universal amplitudes and u is the scaling field 
associated with the leading corrections. With t and B appropriately scaled (so that, 
in particular, the line H c2 (T) becomes t + B — 0), the scaling function C would be 
universal. 

This scaling form has some familiar consequences. In the limit B —> with t > 0, 
we can use (4.6) to find 

C»ci- c 2 t~ a , (5.8) 

with appropriate (non-universal) constants C\ and c 2 . Since a is negative, the spe- 
cific heat rises to a cusp at the critical point t — 0. Below T c in zero field, a 
real superconductor is presumably in its Meissner phase, with (0) 7^ 0, where our 
approximations are not valid. The field dependence at T = T c or t = is given by 

C = d - C 2 B- a,2v C{B 1 - ll2v , uB ul2 ) . (5.9) 

For small fields, we have the asymptotic power law behaviour C ~ c\ — c' 2 B~ a ^ 2u . 
Corrections to this asymptotic behaviour now involve both B 1 - 1 ^ and B"l 2 . In 
the Hartree approximation, it happens that 1 — l/2v = lu/2 = (4 — d)/2 (=1/2 
in three dimensions). A real superconductor, however, is probably characterized 
by the exponents of the three-dimensional XY model, for which w/2 » 0.4 and 
1 — l/2v 0.25, and the leading correction would be that involving B 1 ~ 1 ^ 2v '. 

Of rather greater interest is the behaviour near the line t + eB = 0, corresponding 
to the mean-field H c2 {T) line. According to our earlier discussion, one should even- 
tually pass from the critical region into a region where the lowest Landau level ap- 
proximation becomes good. Within the Hartree approximation, the two-parameter 
scaling form (5.4) is exact, and reduces to the LLL scaling form (5.3) in the limit that 
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5 ~ B w l 2 is large. In general, corrections to the asymptotic, one-parameter critical 
scaling may have many contributions beyond those associated with the scaling field 
u indicated in (5.7), and all of these might be important in the the region between 
the critical and LLL regimes. It is interesting to speculate, on the other hand, that 
this may not be so, and that a crossover to LLL behaviour can be described by the 
two-parameter scaling function C(9,5). This requires that, when 5 is large, 

C(9, 5) « 6 0,/u ''C Ij ll(96~ {2v *~ 1)/uv ) , (5.10) 

where the argument 9S~^ 2v<t '~ 1 ^ u ' v coincides, up to a constant, with the argument 
y — (t + eB) I '(Xf eB)^ appearing in (5.3). We do, of course, find such behaviour in 
the Hartree approximation. 

Unfortunately, mapping out a two-parameter scaling function experimentally to 
test the behaviour suggested in (5.10) would be extremely difficult. Theoretically, 
it is important to note that the crossover mechanism exemplified by (5.10) is quite 
different from that expected at a multicritical point. In the mean field theory of type 
II superconductors, the amplitude of the order parameter in the Abrikosov vortex 
lattice vanishes continuously at the line H c2 (T), and this is sometimes described as 
a second-order phase transition. It is possible to speculate that the specific heat, for 
example, should exhibit an anomaly along this line, perhaps governed by a critical 
exponent a'. In that case, the transition at (T, B) = (T c , 0) would be a multicritical 
point, and one would expect the one-parameter scaling function to incorporate the 
anomaly: 

C(0,O) ~ W a ' (5-11) 

as 9 — > 0. If such an anomaly exists, it is too weak to be resolved by any experiment 
known to us. 
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6. Scaling of the conductivity 

Electrical transport properties of a superconductor can be investigated by using a 
time-dependent Ginzburg-Landau equation to describe the dynamics. A method of 
calculation is described in detail by Ullah and Dorsey 20 who use a Hartree approx- 
imation to study scaling behaviour in the LLL regime. Using essentially the same 
method, we have calculated the conductivities a\\ and a± corresponding to a current 
parallel or transverse to the applied magnetic field. We find that both conductivities 
can be written in the two-parameter scaling form 

a a = B-V»'- d >' 2 S a {e,6) , (6.1) 

with a dynamical exponent z = 2, but with different scaling functions «S|| and <Sj_. In 
the limit of large 8, the LLL scaling properties of the two conductivities are different, 
however, and we find 

« Slll,\\(v) , (6.2) 
where S LL l,\\ = const x T LLL (y)~ l/<t> , while 

a ± » B-^ 2 - d ^ 2 S LLLt± (y) , (6.3) 

with S LLLt ± = const xr LLL (y)~ ( - z+2 ^ d ^ 2 . The results (6.2) and (6.3) agree with those 
quoted by Ullah and Dorsey for d = 3 and d = 2, except that their 2-dimensional 
result for cry is quite different from (6.2), having the form cry ~ B~^ 2 SLLL,\\{y)- 
We are unable to account for this discrepancy in detail. No doubt, however, it has 
to do with the fact that, whereas we have considered an isotropic rf-dimensional 
material, Ullah and Dorsey deal with a layered 3-dimensional system, and obtain a 
two-dimensional limit by taking a large interlayer spacing. 
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7. Comparison with experiment 



While the Hartree approximation may well provide a useful guide to the scaling 
behaviour to be expected in real high-T c materials, it is numerically rather inaccu- 
rate, having critical exponents a — — 1 and v — 1 in three dimensions, compared 
with a ~ —0.01 and v pa 0.67 for the 3-dimensional XY model, which might be 
expected to characterize real materials. Similarly, we do not expect the scaling 
functions to be numerically accurate, and do not present detailed computations. 

Estimates of some gross features of the phase diagram from the Hartree approx- 
imation may, however, have some significance. To obtain such estimates, we need 
to identify the parameters of the model in terms of measurable quantities. We use 
the customary means of identifying the parameters in (2.1), which is explained, for 
example, by Tinkham 21 . These must be treated with caution, however, since they 
are based on the assumption that the entire phase diagram of the superconductor is 
described by mean field theory, which is not actually the case. In SI units, we find 
that the constraint equation (2.8) in 3 dimensions takes the form 

where $ = 2.07 x 10~ 15 Wb is the flux quantum, // is the magnetic permeability, 
which we take to be that of free space (/x = 47r x 10~ 7 Hm _1 ) and k is the Ginzburg- 
Landau parameter. The quantity £ is a characteristic length, of the order of the 
zero-temperature coherence length. 

From the first two terms on the right hand side of (7.1), we find the slope of the 
HdiiT) line as 
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where £o is measured in A and T c in K. We can also obtain an estimate for the field 
strength, denoted by B LLL in figure 1, at which the line H c2 (T) enters the region of 
LLL scaling. At this point, the first two terms of (7.1) cancel, and the argument of 
/ is equal to a value e, for which /(e) ~ ^Jixe^ 1 ! 2 . We thus find 

BlLL ~ ($ e) 3 10 ^ ' (7 ' 3) 

where Bill is measured in tesla. There is no clearly defined value of e which will 
guarantee that the LLL approximation is good. Certainly, e must be smaller than 
the value of approximately 0.3027 for which /(e) = 0. Numerically, we find that f(z) 
can be reasonably well approximated by a function of the form f(z) = f§{z)z~ x l 2 for 
< z < 0.1, where fo(z) is an amplitude which varies from /(0) = y/ir to /(0.1) = 1, 
so we surmise that LLL scaling might be good for z < e — 0.1. 

Detailed measurements of the fluctuation specific heat of YBaCuO have recently 
been reported by Overendet at 22 . For this material, T c = 92K, and we may reason- 
ably take £ ~ 10 A and k ps 100. We then estimate the slope of the H c2 (T) line 
as dB/dT ps —4 TK _1 and, using values of e between 0.1 and 0.2, expect that LLL 
scaling might set in with an applied field B LLL between 10 and 100 T. In fact, the 
measurements of Overend et al indicate that LLL scaling fails at all fields up to 8 T. 
Junod et al 6 have recently reported measurements on the specific heat of YBaCuO 
up to 20 T. For fields greater than 2.5 T, they find that the LLL scaling form col- 
lapses their data onto a common curve, but only in a narrow range of temperature 
near T c2 (B). In the neighbourhood of the peak a little below this temperature, 
their data does not scale at all, in marked contrast to conventional superconductors, 
where the LLL scaling region extends well below this peak 13 ' 16 . For this reason, we 
do not think that their data is really consistent with LLL scaling. There is, however, 
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some indication that the scaling improves at fields of the order of 20 T. Similarly, 
earlier specific heat measurements by Welp et al 17 and by Inderhees et al 4 are at 
best consistent with LLL scaling only in a very narrow range of temperature. 

For a conventional material with, say, T c = 10 K and k — 10 our estimate 
of Bm would be of the order of 10~ 4 T. Thus, the schematic phase diagram of 
figure 1, together with the estimate (7.3) of B LLL is in reasonable accord with the 
observation that LLL scaling appears to work well for conventional superconductors 
over a wide range of fields, but seems not to work for high-T c materials except, 
perhaps, in very large fields. 

Following earlier authors 2 ' 4 , Overend et al attempt to fit their specific heat data 
to a critical-point scaling expression of the form 

C = Ci - C 2 B~^C (^P) , (7.4) 

using the exponents of the 3-dimensional XY model, and find excellent agreement for 
fields up to 8 T, apart from some rounding in very low fields, which is attributable to 
finite-size effects. They find, however, that equally good agreement can be obtained 
using the scaling variable (T — T c2 (B)) / B l / 2v which appears in (5.7), provided that 
the slope of the upper critical field line —dB c2 /dt is greater than about 5 TK _1 . 
This is consistent with the crude estimate of 4 TK _1 obtained above, and with the 
estimate of about 7 TK -1 obtained by Palstra et al 23 by extrapolating the transport 
entropy of vortex motion to zero. However, it would not be consistent with the slope 
of about 1.8 TK -1 required by Junod et al 6 to optimize their fits to LLL scaling. 

The appearance of the scaling field T — T c2 (B), or of t + eB in (4.2) does not seem 
to be an artifact of the Hartree approximation used here. It arises simply from the 
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fact that the Landau eigenvalues are k 2 + (2n + l)eB + t, so that t always appears 
in the combination t + eB, and this would appear to be true quite generally. Within 
the Hartree approximation, we see from (4.5) that the scaling variable 9 ~ (T — 
T C 2(B)) j 'B x l 2v is itself a function of B/t or of where £b is the field-dependent 
coherence length, rather than the zero-field coherence length which appears in the 
gaussian approximation. The difference between the two variables (T — T c (ti))/B l l 2v 
and (T - T c2 {B))jB x l 2v is proportional to B^ 1 / 2 ^ « B - 25 if v is taken to be the 
exponent of the 3-dimensional XY model. Since this correction varies rather slowly 
with B, the distinction between the two scaling variables is probably difficult to 
discern by optimizing the collapse of data onto a common curve. 

8. Conclusions 

We have used a simple Hartree approximation to investigate the critical-point and 
lowest-Landau-level scaling properties of high-temperature superconductors. Our 
principal conclusions are summarized in figure 1, which indicates that critical-point 
scaling is to be expected in the neighbourhood of the zero-field transition, while 
lowest-Landau-level scaling is restricted to a high-field region near the upper crit- 
ical field line H c2 (T). In the critical region, the appropriate scaling variable is 
(T - T c2 (B))/B l / v rather than the variable (T - T C (Q))/B 1 ' 2u which arises in the 
gaussian approximation. While the Hartree approximation does not yield accurate 
values for critical exponents or scaling functions, it suggests that a minimum applied 
field of between 10 and 100 T is required to observe lowest-Landau-level scaling in 
materials such as YBaCuO. These conclusions appear to be consistent with current 
experimental observations. 
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Figure Caption 

Figure 1. Schematic phase diagram of a high-temperature superconductor in the 
neighbourhood of its critical point. Shaded regions are those in which critical point 
or lowest-Landau-level scaling should be observed. The line t = —eB corresponds 
to the upper critical field, while the line (0) ~ is a notional one, below which the 
approximation (0) xs might be expected to fail. The curve t/2eB = e represents a 
criterion for the validity of the lowest Landau level approximation explained in the 
text. 
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